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A deposition process with particles having realistic intermediate stickiness is studied in 2 + 1 
dimensions. At each stage of the deposition process, for any given configuration, a newly depositing 
particle gives rise to allowed set of configurations that are vastly larger than those for deposition of 
a mixture of purely non-sticky (random like) and purely sticky (ballistic like) particles. We obtain 
scaling behaviour and demonstrate collapse of scaled data for surface width and porosity. Scaling 
of conductivity, when a porous structure thus formed, is saturated with conductive fluid, e.g. brine, 
is studied. The results obtained are in good agreement with Archie’s law for porous sedimentary 
rocks. 


I. INTRODUCTION 

The growth of surfaces in different dimensions and on different substrate geometries, finds applications in several 
areas of science and technology, including physics, chemistry, biology, geology, chemical engineering and material 
science. Though the applications are diverse, physical properties of growth processes, such as, nature of roughness, 
porosity and their dynamic scaling behavior depend on few basic entities, such as, dimensionality, geometry and 
underlying symmetry of the problem. Hence, these systems can be classified into a few universality classes. Extensive 
theoretical and experimental study have been undertaken in these areas. The theoretical study of dynamic scaling 
behaviour of surfaces generally follows two pathways, extensive simulations of discrete models Em, and study of 
relevant stochastic differential equations obtained from phenomenological consideration SI- Another relatively 
recent approach, involves obtaining difference equations from microscopic deposition rules for the discrete models and 
deriving relevant stochastic differential equations by limiting process using various regularization techniques. EM- 

The simplest growth process is studied by considering a single species of particles, either completely non-sticky 
or completely sticky, descending on a one-dimensional or a two dimensional substrate. These are called random 
deposition (RD) and ballistic deposition (BD) respectively. For more realistic systems, one needs to consider the 
possibility that, a single type of particle may have an intermediate stickiness and in one deposition process, several 
such species may be involved. 

In order to study intermediate stickiness, Wang, Cerdiera et.al. have studied models with two types of particles, 
some random-like, and others ballistic-like [lll - [l4| . However neither species are allowed to have intermediate stickiness. 
Horowitz, Albano mm have studied growth models in which each incoming particle may behave either as non-sticky 
with probability p or as completely sticky with probability (1 — p). However, neither mimics the possibility that a 
particle may be partially sticky. Thus at each contact with the surface, it may have a fixed probability of sticking 
(0 < p < 1), and a fixed probability (q = (1 —p) < 1) of continuing its journey till it settles somewhere on the surface. 
Study of such model in (1 + 1) dimension was proposed in an earlier workl_17j. In the present article, we extend the 
study to (2 + 1) dimensions. A related model on a one-dimensional substrate , with next nearest neighbor sticking 
was studied by Banerjee et al. fl8l ]. 

In this physically realistic model, the incoming particle may come in contact with several points on the surface, 
and its final position of deposition forms a vastly larger ensemble than that considered in the former studies (Tll - [l6| . 
In case of (2+l)-dimension, i.e., for growth on a two dimensional substrate, the present model may find application 
in the formation of sedimentary rocks and fabrication of nano materials. 

The structure of pores in depository rocks is of great importance in rock geology (petrology) and oil exploration. 
Surface roughness is measured in terms of the standard deviation, which is the square root of the second central 
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moment of the height distribution and is defined as, 
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where, h(i,j,t) is the height of the (i,j )-th site at any instant t, L is the system size and () is the average. Though 
this is an important quantity of interest, it is far less informative than the distribution of height itself. Knowledge 
of frequency or probability distribution of height is equivalent to the knowledge of all moments [l9|. Similarly, the 
porosity gives an average information about the nature of the pores. It cannot, however, tell us if the pores are clumped 
together, or are more scattered throughout the allowed volume. Moreover, it cannot distinguish between situations 
where pore clusters are elongated longitudinally or transversely. Additional information about conductivity may throw 
some light on the above mentioned geometry, though in a somewhat qualitative manner. For a given porosity, higher 
conductivity implies more clustering of pores, predominantly in the longitudinal direction, rather than in circuitous 
paths having long horizontal parts. In addition, in real life experiments in geology, this conductivity is more easily 
measured than other measures of pore geometry. Thus the study of conductivity, and its scaling with porosity is 
important for better understanding of the porous structure. In this article we obtain surface width, porosity and 
conductivity, and their scaling behavior in (2+l)-dimensional growth on a flat substrate. 

In RD the individual columns grow independently of each other without any bound and thus roughness of the 
interface width never saturates. No voids are present within the composite [lj. Correlations can be introduced by 
making the particles sticky. The model representing such a system is called ballistic deposition (BD) model, where 
an incoming particle sticks to the first point of contact it encounters with the surface while falling down vertically 
towards a randomly chosen site on the substrate 111 The stickiness in BD is extreme, the particle must stick 
at the very first contact and is given no option to slide past the first point of contact. In real systems, one may find 
a particle stickin g to a site after sliding past a few points of contact on its vertical journey. In this article we extend 
our earlier work 17] of such a realistic process, to surfaces growing on a two-dimensional substrate and building a 
three-dimensional structure with voids. This porous structure may fill up with a conducting fluid, such as, brine, as in 
the case of sedimentary rocks. The conductivity of the structure depends on the conductivity of the brine, the relative 
amount of pores, and the geometry of the pores in the 3-d structure. To begin with, we take the specific conductivity 
of brine as a constant. The porosity and geometry of the pores both depend on the stickiness and possibly the size of 
the substrate. The variable stickiness of particles is modeled using a parameter 1 > p > 0. A particle dropped on to 
the substrate sticks to the first surface it encounters with a probability p and continues on its downward journey with 
probability (1 —p). The probability that it will deposit at the very next surface it encounters isp(l— p). It will continue 
to the next lower position with probability (1 — p)(l — p). Thus, if a site is selected with taller nearest neighbor, a 
newly arriving particle can deposit at any one of the successive positions 1,2,3... shown in Fig. [TcJ with probabilities 
p,p( 1 — p),p( 1 — p) 2 ■ ■ ■ respectively. In this model, p = 0 corresponds to RD (Fig. [Tall and p = 1 corresponds to 
BD (Fig. Ilbl) . while 0 < p < 1 represents intermediate stickiness. It may be noted that the present model does not 
allow sticking at corners or on edges, that is it disallows sticking to next nearest neighbors. In the present model 
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FIG. 1: Allowed positions in (a) RD (b) BD (c) GBD 


the surface width depends on the sticking probability p and system size L, in addition to time t, measured in terms 
of the average number of layers deposited during that time. The logarithmic plot of surface width W versus time t 
shows four distinct regions. There is an initial random like growth region (GR-1), followed by a non-KPZ like growth 
(GR-2) and then a KPZ growth (GR-3) followed by an eventual saturation W(L,p,t ) —>• W sat {L,p). Similar feature 
was also reported earlier in one dimension E3, and in systems where particles may stick to next nearest neighbors 
fist . With the introduction of probability of sticking p, the Family-Vicsek scaling relation W(L,t) ~ L a f(t/L z ) is 
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modified and the following dynamic scaling relation obtained for the growth and saturation regions, 


W(L,p,t) ~ L a p~ a ' F 



( 2 ) 


where F(x ) is scaling function that satisfies F( oo) ~ constant and F(x) ~ x& for small x. The scaling exponents are 
determined and an excellent collapse of scaled data is obtained using those exponents. 

We also study the porosity of the bulk of size L 3 just beneath the active region, which is signified by the volume 
inaccessible to the new incoming particles. It is found that the porosity reaches a constant value almost as soon as 
a bulk volume of dimension L 3 forms below the active region. We find a scaling relation between saturated porosity 
p g , system size L and probability of sticking p as, 

Ps ~ p a L b . (3) 

The porous structures thereby formed are further investigated for their conducting properties. The saturated con¬ 
ductivity a s is found to obey a scaling relation 

cr 8 ~ p m L n . (4) 

The corresponding exponents are calculated from our simulational results. In addition we observe that (a s /L) depends 
on p s as (<j s /L) ~ /?{ with / = 2.02, which is in good agreement with Archie’s law, an important empirical law in 
geophysics 121 It23| . 

In our simulations of the present generalized deposition model, two independent random number generators were 
used, one for selecting a site on the growing surface and, another to determine whether a particle will stick at a 
particular location for a chosen value of the sticking probability. These two random number generators are chosen to 
ensure that they are completely independent and uncorrelated. The reliability of the random numbers used is verified 
by the x~test and absence of repetitive subsequences or looping for the maximal set of random numbers drawn for 
both the sets. 


II. GENERALIZED BALLISTIC DEPOSITION (GBD) - VARIABLE STICKINESS 

The present model is a modification of ballistic deposition to represent realistic sticky particles. The model is 
studied in (2+l)-dimensions. A particle is allowed to descend vertically towards a randomly chosen site on a two 
dimensional substrate. If the selected site is higher than its nearest neighbors, the particle simply deposits on top of 
the column at that site. However, if the chosen site has a taller column of particles as its nearest neighbor, then the 
new particle sticks to the first occupied site it encounters if the value of p is larger than a random number generated 
from a uniform distribution between 0 and 1. Otherwise, it slides down vertically to the next occupied site with 
probability (1 — p). At this site the particle may stick with probability p( 1 — p) or continue its further descent with 
probability (1 — p) 2 , and so on, till it reaches the bottom. Thus if the chosen site has a nearest neighbor with column 
height taller by n layers relative to it, the probabilities of the arriving particle sticking to the successive particles of 
the nearest neighbor column from top are given by, 

P(l)=p, P(2) = p(l — p),... P(k)=p(l-p^ k ~ 1 \ (5) 

The probability that the particle slides past the preceding (n — 1) occupied neighbors, and lands at the lowest possible 
position is given by, 

n —1 

P(n) = l-^P(fc) = (l-rt (n - 1) . 

fc=l 

This describes a proper stochastic process. The total probability of a descending particle sticking to one of the allowed 
position is Y^k= l -P(&) = 1- It must be noted that sticking to corners or edges are not allowed. Only surface sticking 
is allowed in this particular model. 


III. RESULTS AND DISCUSSIONS 
A. Scaling of surface width 

Simulations have been performed for several system sizes and various probabilities of sticking. We present 
the analysis of data for system sizes L = 64,128,256,512 and values of probability of sticking p = 
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1,0.8,0.7,0.5,0.25,0.125,0.0625. The p = 0 gives us the random limit and p = 1 is the ballistic deposition. The 
logarithmic plot for surface width and time shows four distinct regions with varying slopes as shown in Fig. [2] It is 
interesting to note that this feature, namely the existence of four characteristic regions, is observed whenever stickiness 
is present, i.e. for all non-zero probability of sticking p, however small. 

The dependence of surface width W on t in log-log scale, in the early submonolayer region (t -C 1) is linear with 
slope 1/2 as in random deposition (growth region 1, GR-1). At later stages of submonolayer growth (growth region 
2, GR-2), t ~ 1”, the surface width shows a steep increase which continues for the first few layers (1 — e < t < 3, 
1 e > 0). With deposition of further layers, the rate of increase in width slows down (growth region 3, GR-3). After 
deposition of a large number of layers, the ensemble average of the surface width saturates. Three different crossover 
times are of relevance. The first crossover time t r corresponds to the change from random growth to region with slope 
greater than 1/2. The second crossover time tk corresponds to time beyond few layers where the slope decreases and 
changes from GR-2 to GR-3. The third crossover time t sat corresponds to beginning of saturation region. 




FIG. 2: Logarithmic plot of interface width and time for 


FIG. 3: Characteristic plot for lnIT sa t vs In t with 


different p for L = 256. 


p = 0.5 for different system sizes. 




FIG. 4: Deviation from the random like behavior at later 
stages of submonolayer regime for L = 256. 


FIG. 5: Dependence of width in later stages of 
deposition with saturation for L = 256. 


The appearance of different growth regions in the present model may be understood as follows. In our model, 
we start from a flat substrate, thus initially almost no two adjacent sites are occupied, hence there is no correlation 
among neighboring columns. Thus for system of all sizes, at the very beginning, when t <C 1, the growth is random 
like, irrespective of whether we allow sticking or not (Fig. [2] and Fig. [3]). The deviation from random like behavior 
begins near the first monolayer, and may continue for few monolayers of deposition. In this region the surface width 
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grows at a rate faster than that in the case of random deposition. As the number of particles deposited at this stage is 
nearly L, due to fluctuation, some short multi-layer columns begin to form. This brings in non-trivial correlations in 
the system, due to possibility of the descending particles encountering occupied neighbors before reaching the bottom 
of their own columns. At this stage, the average height of the surface is small, and even a few particles sticking to a 
higher location instead of reaching the bottom of a column, makes a significant relative change in width. Thus the 
rate of growth of surface width in this region is higher than that for RD (Fig. 0J. 

With further layers of particles being deposited, we reach a third growth region GR-3, where the rate of increase 
in surface width slows down (Fig. [5]). The average height and the interface width are large in this region. The deep 
crevices are encountered by descending particles, and if they are sticky, they can stick to a side wall, thus filling up 
the crevices much faster than for RD, where the crevices need to be filled from bottom up. At an even later time, the 
above mentioned smoothening effect starts dominating, and the surface width finally saturates. The saturated width 
depends both on the system size L and sticking probability p. 

For a given value of p , the saturated width W sa t and t sat increase with system size L and for fixed L, they decrease 
with increase in probability of sticking p. This decrease is more pronounced for lower values of stickiness parameter, 
i.e., p < 0.5. The KPZ-like growth region GR-3, and the saturation region follow scaling relation stated in Eqj2] 
for which the exponents can be evaluated. With the increase in the probability of sticking p , the saturation is at 
lower values of interface width. This dependence is found to be of the form p~ a with a' = 0.20846. Further, this 




partially scaled width lnW p a depends on the system size as L a with a = 0.3224 (Fig. [6]). The plot of scaled width 
(In Wp a / L a ) versus Inf is shown in Fig. 0for p = 0.8,0.5,0.25,0.125 for system sizes L = 64,128,256,512. The 
figure shows collapse of scaled data in the saturation region. We observe that the KPZ-like growth region is more 
pronounced, and occurs over a wider time interval for larger system sizes, as the correlation effects take a longer 
time to bring in the saturation. However, with increasing probability of sticking p , the correlation effects are more 
dominant and hence saturation kicks in much earlier, thereby resulting in a shorter KPZ-like growth region. The 
scaling for the growth region with respect to sticking probability is shown in Fig. [51 and the complete scaling is shown 
in Fig. [9l using the our calculated values for exponents z = 1.72348 and z' = 1.069246. The slope of the fully scaled 
width versus time, in the growth region GR-3 is determined as /3 = .19. Thus in growth region we get, 


W(L,p, t) ~ L a p- a F 
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Whence we identify the following relation connecting the obtained scaling exponents, 
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where a = .3224 ~ 1/3, 0 = 1.7234 ~ 5/3, a' = .20846 ~ 1/5, z' = 1.069 ~ 1 and /? = .19 ~ 1/5. 


( 6 ) 

( 7 ) 
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B. Porosity and Conductivity 

In this section, we study the dependence of porosity on system size L. Porosity p is defined as the number of vacant 
sites within a cubic volume of side L , just beneath the active surface of the deposit. For (2 + l)-dimensional systems, 
it is found that the porosity reaches a saturation value by the time a volume L 3 is formed below the active layer. 
Thus, there is no measurable variation of the porosity with time. The dependence of the porosity on the system size 
L is also not very significant (Fig. [Hi]). The deposit becomes more porous for higher sticking probability p^> 0. The 
saturated porosity p s scales with p and L as, p s ~ ( p a L b ) with a = 0.182920 and b = 0.00478. 
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FIG. 10: Variation of p s with system size L in log-log 
scale for different values of p. 


FIG. 11: In p s /L b versus lnp. 


The conductivity a depends on the porosity p as well as the geometry of the pore structure. In order to study 
the dependence of a on p , an inactive sample of dimensions L 3 is considered where no more deposition can take 
place. The solid deposited particles are insulating whereas the voids are considered as conducting when filled with a 
conducting liquid, such as brine. Thus, the connected pore space is first identified using Hoshen-Kopelman algorithm 
MM- Laplace’s equation is solved for the potential distribution to obtain the conductivity. The discrete version of 
the Laplace’s equation V 2 V ( x , y,z) = 0 is given by 

V— i,j,k T Vi+i,j,k , Vp —qfc 2Vj, i ?,fc V , Vi,j,k— l V n 

(Air) 2 + (A^p + (Ai) 2 


( 8 ) 
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FIG. 12: Variation of p s with p in log-log scale for 
different values of L. 


FIG. 13: In p s /p a L b as a function of InL and lnp. 


We consider (Ax) = (A y) = (A z) = 1, as deposition consists of unit cubes. Thus, Eq [8] becomes, 


v. (”ti) = a w (n> -t -v yn> +vy i i 4-v v,l) +vy i> +vy i> , i 

v i, j,k at v i— l,j,k ' v i,j—l,k ' v i,j,k— 1 ' v i-\-l,j,k ' v i,j-\-l,k ' v i,j,k +lJ 
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where, N is the number of nearest neighbor vacant sites. The boundary conditions are V(z — 0) = 0, V{z = L) = 1. 
The initial condition obtained by setting values of V^ for n = 0. We have started with V^ = ( k/L ) for all vacant 
sites at a height k. Numerically, the steady state is said to be achieved if (| 1 [A ,l+1 ) — < e), where e is the 

required accuracy. If the steady state condition is achieved after r iterations then V r is the steady state potential. 
The conductivity is also found to be a constant cr s for a fixed system size L and p. For a given L and p , a reaches a 
steady state value a s almost simultaneously with porosity p. a s depends on system size L and sticking probability p 
as <T s ~ p m L n with m = .368496 and n = 1.02323 (Fig. [fdland Fig. fl5l) . Finally, the dependence of a s on p s is found 
to satisfy Archie’s law of the form a s ~ p{ with / = 2.02 as shown in Fig. [181 An excellent collapse is obtained for 
a s as shown in Fig. flUland Fig. [T7]with m = 0.368496 and n = 1.02323. 
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FIG. 14: Variation of a s with system size L in log-log 
scale for different values of p. 


FIG. 15: Variation of cr s with p in log-log scale for 
different values of L. 
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FIG. 17: \na s /p m L n as a function of InL and lnp 



FIG. 18: Variation of a s /L on p 


IV. CONCLUSION 

We have studied a generalized ballistic deposition (GBD) model with deposition of particles having intermediate 
stickiness in (2+1) dimension. Surface width, porosity and conductivity, when such structure is saturated with 
conducting liquid, show scaling behaviour with both system size L and sticking probability p. Scaling relation of 
surface width is studied and correct scaling exponents are determined. The deposition process leads to porous 
structures. Scaling of saturated porosity and conductivity with system size and sticking probability are also studied 
and the corresponding scaling exponents are calculated. The scaling exponents thus obtained agree with Archie’s law 
for the dependence of the conductivity on porosity. A study of generalised deposition with more species of particles 
having varying intermediate stickiness 1 > pj. > 0, with k = 1,2,3, ...zz, where v is the number of such species, may 
be of interest in depository rocks. Such a study is in progress and the results will be reported elsewhere. 


[1] A. L. Barabasi and H. E. Stanley, Fractal Concepts in Surface Growth (Cambridge University Press, Cambridge, 1995). 

[2] P. Meakin, Physics Reports 235, 189 (1993). 

[3] F. Family and T. Vicsek, J. Phys. A: Math. Gen. 18, L75 (1985). 

[4] F. Family, J. Phys. A: Math. Gen. 19, L441 (1986). 

[5] S. F. Edwards and D. R. Wilkinson, Proc. R. Soc. Lond. A 381, 17 (1982). 










9 


[6] M. Kardar, G. Parisi, and Y.-C. Zhang, Phys. Rev. Lett. 56, 889 (1986). 

[7] D. D. Vvedensky, Phys. Rev. E 67, 025102 (2003). 

[8] D. D. Vvedensky, Phys. Rev. E 68, 010601 (2003). 

[9] C. A. Haselwandter and D. D. Vvedensky, Phys. Rev. Lett. 98, 046102 (2007). 

[10] C. A. Haselwandter and D. D. Vvedensky, Phys. Rev. E 76, 041115 (2007). 

[11] W. Wang and H. A. Cerdeira, Phys. Rev. E 47, 3357 (1993). 

[12] W. Wang and H. A. Cerdeira, Phys. Rev. E 52, 6308 (1995). 

[13] H. F. El-Nashar and H. A. Cerdeira, Phys. Rev. E 61, 6149 (2000). 

[14] H. F. El-Nashar, W. Wang, and H. A. Cerdeira, J. Phys.: Condens. Matter 8, 3271 (1996). 

[15] C. M. Horowitz and E. V. Albano, Phys. Rev. E 73, 031111 (2006). 

[16] C. M. Horowitz and E. V. Albano, J. Phys. A: Math. Gen. 34, 357 (2001). 

[17] B. Mai, S. Ray, and J. Shamanna, Arxiv: 1410.6600 (2014). 

[18] K. Banerjee, J. Shamanna, and S. Ray, Phys. Rev. E 90, 022111 (2014). 

[19] B. Mai, S. Ray, and J. Shamanna, Eur. Phys. J. B 82, 341 (2011). 

[20] F. Family, Physica A: Statistical Mechanics and its Applications 168, 561 (1990). 

[21] G. Archie, American Association of Petroleum Geologists Bulletin 31, 350 (1947). 

[22] G. Archie, American Association of Petroleum Geologists Bulletin 34. 

[23] K. Verwer, G. Eberli, and R. Weger, American Association of Petroleum Geologists Bulletin 95, 175 (2011). 

[24] J. Hoshen and R. Kopelman, Phys. Rev. B 14, 3428 (1976). 

[25] D. Stauffer and A. Aharony, Introduction to Percolation Theory (Taylor and Francis, London, 1992). 



